


> 



X 



DESY 96-096 |hep-ph/9712546 



WUE-ITP-96-39 December 1997 



The Evolution of Unpolarized Singlet 
bo Structure Functions at Small x 



Johannes Bliimlein" and Andreas Vogt^ 



^ , ''DESY-Zeuthen 

<3\ ■ Platanenallee 6, D-15735 Zeuthen, Germany 



Ph, ^Institut fiir Theoretische Physik, Universitdt Wurzburg 



Oh! Am Huhland, D-97074 Wurzburg, Germany 



Abstract 

A systematic study is performed of the impact of the various resummed small- a; contribu- 
tions to the anomalous dimensions and coefficient functions on the evolution of unpolarized 
structure functions in deep-inelastic scattering. The proton structure functions F2 and Fj^ 
as well as the photon structure function F2 are considered together with the correspond- 
ing parton densities. The general analytic solution of the evolution equations in Mellin-A^ 
space is derived, and different approximate solutions are compared. Potential effects of less 
singular small-x terms in the anomalous dimension and coefficient functions are discussed. 



PACS numbers : 13.60.Hb, 12.38. Cy, 12.38bx 



1 Introduction 

One of the central questions in the theory of deep-inelastic scattering (DIS) structure functions is 
that of their behavior at small values of Bjorken-x. The HERA experiments [jl]] have performed 
detailed measurements of the structure function -F2(x, Q^) down to values of x ^ 10~^ and have 
presented first results on Fl{x,Q'^) in the range x>10~^ 0. F2 rises even at small values of 
Q^ ~ 2 GeV^ approximately as x"*^'^. At low scales Q^ the behavior of structure functions 
cannot be dealt with by means of perturbative QCD due to the size of the strong coupling 
constant ^^(Q^)- For large virtuahties, on the other hand, a perturbative description of the 
scaling violations is possible if a factorization can be achieved between the non-perturbative 
input distributions and the evolution kernels which can be calculated perturbatively. 

Throughout the present paper we will consider deep-inelastic scattering in the range where 
it is dominated by the light-cone singularities in the Bjorken limit, Q^, s -^ 00, x = Q'^/{sy) = 
const., with Q^ the 4- momentum transfer squared, s the center-of-mass energy, y = 2P-q/s. The 
ultraviolet singularities of the operators emerging in the light-cone expansion Q] are associated to 
renormalization group equations (RGE) which describe the evolution of the structure functions. 
Under these conditions the evolution kernels are given by the anomalous dimensions of the light- 
cone operators. For the leading-twist contributions considered in the following, the expectation 
values of the operators are related to the parton distribution functions. This picture holds, in 
principle, down to the region of small values of x. The anomalous dimensions and coefficient 
functions, however, may receive large low-a; contributions of the type 0] 

a^f^^X ^ l^.,.ln^::(lM, 

'\N-lJ X ' (k-iy. ' 

where N denotes the index of the Mellin-transformation 

M\f{x)]{N) = f dxx^fix) . 
Jo 

The x-range in which this representation is applicable can only be determined by explicit calcu- 
lations 1^ and depends on both the behavior of the partonic input densities at a starting scale 
Qq and on the structure of the anomalous dimensions and coefficient functions. Note that the 
present picture does not necessarily yield a description of the structure functions as well in the 
Regge limit s -^ 00, Q"^ = const., since both limits lead in general to different results. This 
can be seen performing these limits in the Jost-Lehmann-Dyson representation |j^ of the 
structure functions. 

In the Bjorken limit the ultraviolet and collinear divergences emerging in the calculation of 
the higher-order corrections can be dealt with applying the corresponding RGE-operators, which 
imply the evolution of the parton densities and the running of the strong coupling constant. 
The resummed small-x corrections in leading |^ and next-to-leading order [8-10] can thus be 
accounted for in a natural way. Since the impact of the resulting all-order anomalous dimensions 
on the behavior of the DIS structure functions at small x does as well depend on the non- 
perturbative input parton densities at an initial scale Qq, the perturbative resummation effects 
can only be studied via the evolution over some range in Q^. This evolution probes also the 
anomalous dimensions and coefficient functions at medium and large values of x due to the 
Mellin convolution between the evolution kernels and the parton densities, cf. ref. [jrT|. Hence 



the small-x dominance of the leading terms over less singular contributions as x — *> in the 
anomalous dimensions and coefficient functions does not necessarily imply the same effect for 
the observables, such as the structure functions. 



In the present paper a systematic study is performed on the impact of the different resummed 
small-x contributions to the anomalous dimensions and coefficient functions on the evolution of 
the singlet contributions to the nucleon structure functions F2{x,Q'^) and Fl{x,Q'^). A brief 
summary of some of our numerical results on the effects of the resummed gluon anomalous 
dimension on the nucleon structure functions has already appeared in [|I^. We extend the 
analysis to the photon structure function F2{x,Q^), for which we derive corresponding results 
in the DIS scheme. 

The paper is organized as follows. The general framework for the evolution of parton densities 
of the nucleon and the photon is outlined in Section 2. In Section 3 the presently known small- 
X resummed anomalous dimensions and coefficient functions [4,8-10] are summarized and the 
numerical coefficients for their expansions in a^ are compiled for the subsequent analysis. The 
issue of less singular terms is discussed in Section 4 guided by the known 2-loop results. In 
Section 5 different methods are derived for the solution of the evolution equations in the presence 
of all-order resummations for the small-x contributions. The numerical implications of the small- 
X resummations on the evolution of the parton densities and the structure functions Ff, F^ and 
F2 are worked out in detail in Section 6. Section 7 contains our conclusions. 

2 The evolution equations 

The twist-2 contributions to the structure functions in inclusive deep-inelastic scattering can be 
described in terms of the QCD-improved parton model. Their scaling violations are governed 
by renormalization group equations which can be formulated to all orders in the strong coupling 
constant. In this section we briefly recall this general framework, which allows for a consistent 
introduction of the small-x resummations into the structure function evolution, and its specific 
application to hadronic and photonic parton distributions. 

2.1 The renormalization group equations 

Among the singularities emerging in the calculation of QCD radiative corrections, only the 
ultraviolet divergences and the initial-state mass singularities require a special treatment in 
inclusive deep- inelastic scattering [|. The former are eliminated by the renormalization of the 
strong coupling constant. The remaining mass singularities, originating in collinear emissions of 
massless partons off massless partons, are removed by mass factorization, i.e., these contributions 
are absorbed into the bare parton densities. For this procedure the structure functions Fi{x, Q^) 
are first written as 

F,(x, Q2) = F. , (^^, asiR''), %, ^, e) ® A(x) . (2.1) 

Here ^^(i?^) denotes the strong coupling constant, renormalized at the scale R in some renor- 
malization scheme, n is an arbitrary mass scale, fk represents the bare momentum distribution 
of the parton species k, and ® stands for the Mellin convolution in the first variable, 

A{x,jjI) 1^ B{x,jjI) = / dxi / dx2 6(x - X1X2) A(xi, fil)B(x2, i4,) ■ (2.2) 

JO Jo 



^The infrared divergences cancel between the virtual and Bremsstrahlung contributions according to the 
Bloch-Nordsieck theorem [Q. The final-state mass singularities vanish due to the Kinoshita-Lee-Nauenberg 
theorem 1l4], as all degenerate final states are summed over. 



Finally e marks the initial-state mass singularities entering the bare partonic structure func- 
tions Fi^k- These functions are then separated into the coefficients functions Cjj and the transi- 
tion functions Fj^, which contain the 1/e pole terms, according to 

Fi^k(x,as{R^),^;,^,e] = C^^^f x,a,(i?^), — , — j ® F^^f x,a,(i?^), ^, — ,ej . (2.3) 

The additional parameter M is the factorization scale. This separation is not unique beyond 
the leading order (LO) of the perturbative expansion, hence T^j and Ci^k are also factorization 
scheme dependent. Combining Eqs. (p. ID and (p.3|), the structure functions Fi{x, Q"^) finally read 



F,{x, Q') = Q, [x, a,{R% ^, |^) ® /, [x, a,{R% ^, |^) (2.4) 

in terms of the renormalized parton densities fj given by 

/j|x,«,(i?2),^,— j =T^k[x,as{R^),^^,^,e^®fk{x). (2.5) 

Two arbitrary scales, R and M, are thus introduced by the renormalization and mass fac- 
torization procedures. These scales are not physical. Hence observables, such as the structure 
functions Fi{x,Q'^) in Eq. ( p.4|) , do not depend on them. Since the ultraviolet and mass singu- 
larities are not related, the conditions 

R'^F^(^^Q') = 0' (2-6) 

M'^^F^i^^Q') = (2.7) 

hold separately and imply two independent renormalization group equations. The first of these 
equations leads to the scale dependence of the running coupling constant, 

"^""^ /?(«.) = -f:«^'A, (2.8) 



d In i?2 ^^g 



where the abbreviation Og = Q;s(i?^) /(47r) has been introduced for convenience. At next-to- 



leading order (NLO) only the first two, scheme independent terms [|l^, [T^ of P{as) are kept 
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A = -CI--CaTf-4CfTf. (2.9) 

The QCD color factors are Cp = {N^ - l)/(2A/'c) = 4/3, Ca = N^ = 3, Tr = 1/2, and 
Tp = NfTpi^, with Nf denoting the number of light quark fiavors. To this approximation, the 
solution of Eq. ( |2.8| ) can be expressed in terms of the QCD scale parameter A^r^ by 

In turns out, moreover, that higher coefficients /3;>2 - /32 and jS^ have been calculated in the MS 



scheme |T^, [T^ - are not required in connection with the presently available small-x resumma- 
tions, see Sect. 5.3. We will therefore employ the relation p.lOp for all our numerical calculations 
in Sect. 6. 



2.2 Hadronic and photonic parton densities 

The second renormalization group equation ( |2.7D leads to the scale evolution of the renormalized 
parton densities fj{x,M'^). Considering first the hadronic case, the relevant parton species are 
the quarks and antiquarks, qj and q^, and the gluon g. It is convenient to introduce flavor 
non-singlet combinations of the quark densities, 

qf = qj±qj~jrT.[ir±Qr], (2-ii) 

Nf 

q""' = Y.[^r-qr]. (2.12) 

and the singlet quark / gluon vector 

9=f^l , ^^i:\qr + qr]. (2.13) 



\ ^ / r=l 

This decomposition decouples the 2Nf + 1 evolution equations as far as possible by symmetry 
considerations alone. For simplicity, we will choose the renormalization and factorization scales 
as B? = M^ = Q^ from now on^ The hadronic evolution equations can then be written as 



dqf{x,Q^ 



' ' ' - P^{x,as)^qf{x,Q^), 

P(x,a,) ® q(x,Q'). (2.14) 



dqix,Q') 



91nQ2 

The evolution of g™' is identical to that of q~ up to NLO. As our present analysis is confined to 
that approximation in the non-singlet sector 0, the corresponding equation has been suppressed 



in Eqs. ( 2.14 ). In general, the splitting functions P , P are given by the infinite series 



oo 



1=0 

The expansion coefficients Pf{x) and p™P° i^x) are, in sensible factorization schemes, subject 
to the sum rules 



dxPr{x) = Q, (2.16) 

fdxxY.P-kt-(^) = ^^ (2-17) 



^° J 



^See refs. [Q Q for studies of the uncertainties in NLO analyses due to the variation of R and M. 

^The generating relations for the resummation of the anomalous dimensions |^ of the ± non-singlet combi- 
nations were derived in ref. |2^. Here the leading small- a; terms are of 0[(as In a;)"]. The effect of these terms 
has turned out to be on the 1% level down to very small values of x |^, g3|. As the non-singlet contributions 
are furthermore suppressed compared to the singlet ones at low x, these resummations are not included in the 
present treatment. 



which are due to fermion number and energy-momentum conservation, respectively. By now all 
unpolarized and polarized entries in Eqs. ( |2.15| ) are completely known at NLO, / = 1. The full 
expressions for their x-dependences can be found in refs. [24-28]. Beyond this order a series of 
integer Mellin moments of P2~{x) and P2{x) has been calculated so far [^ . 

We now turn to the parton densities of the real photon. The photon is a genuine elementary 
particle, unlike the hadrons. Hence it can directly take part in hard scattering processes, in 
addition to its quark and gluon distributions arising from quantum fluctuations, q^i^x^Q"^) and 
g'^{x,Q'^). Denoting the corresponding photon distribution in the photon by V^lXjQ"^), the 
evolution equations for these parton densities are generally given by [EO] 



dg^ 



^ fc=i 



N 



dg 



7 



ainQ2 

aing2 



aemPg-y ® T ^ + a J 2 ^ Pg,, ® q^ + Pgg ^gH, (2.18) 

*- fc=i ^ 

Nf 

aem^77 ® T"^ + Oem] 2 ^ P^g^ (g) q^ + P^^g ^g^'i . 
^ fc=l '' 



Here Oem = a/(47r) with the electromagnetic coupling constant a ~ 1/137. The antiquark 
distributions do not occur separately in Eqs. (|2.18|) , as q^{x,Q'^) = q^{x,Q^) due to charge 
conjugation invariance. The generalized splitting functions read 

rij[X, Clem} Ois) ^ /-^ ^enfls ^ij V-^J 5 \Z,.l\)) 

l,m=0 

with Pq^Qf. being the average of the quark-quark and antiquark-quark splitting functions. 

Usually calculations involving the photon's parton structure are restricted to the first order 
in a <^ 1. In this approximation all m 7^ terms in Eq. (|2.19|) can be neglected, since q^ and g"^ 
are already of order a. This reduces the functions P^ to the usual QCD quantities Pij{x,as), 
and P^q^ and P^g drop out completely. Moreover one has P^^ oc S{1 — x) to all orders in as, as 
real photon radiation from photons starts at order a^ only. Thus the last line of Eq. (|2.18|) can 
be integrated immediately, at leading order (LO, / = 0), for example, resulting in 

Tl^oI^, Q') = '^(1 - a;) [1 - 4a, ( ^ ej^ In ^ + const.)] , (2.20) 



fc=i 



Ql 



where Cq^, represents the quark charges, and Qq is some reference scale for the evolution. Only 
the 0(1) part of V^ affects the quark and gluon densities at order a, as well as any observable 
involving hadronic final states like F2 ■ Therefore, after decomposing into the singlet and non- 
singlet parts as before, one obtains the inhomogeneous evolution equations 



'r+f^ n2^ 



dq]^{x,Q 



Q^^g2 = k+{x,as) + P+{x,as)^q]^{x,Q^), 
■)2^ 

-l2^ 



dq^ix,Q' 



d\nQ^ = H^^(^s) + P{x,as) ® qHx,Q') . (2.21) 

Here we have switched to the conventional notation for the photon-parton splitting functions 

00 
k^{x,as) = ^aeniai/^i^(a;) , 

1=0 



k{x,as) = (^'&'li]=f:aerna[ki{x). (2.22) 

These splitting functions are presently also known up to NLO (/ = 1), see refs. [30-32]. 

2.3 Factorization scheme transformations 

As stated above, the separation between the coefficient functions and the splitting functions is 
not unique beyond the leading order. In this subsection, we derive the general factorization- 
scheme invariants and specify the schemes for our subsequent numerical calculation. We study 
the photonic case, as the hadronic problem forms a subset hereof. One way to carry out this task 
is to introduce a fictitious second, 'gluonic' structure function, where the gluon density enters at 
order a°, but the quarks only at order a], opposite to the situation with the real electromagnetic 
current, see refs. [^, |3^. This can be achieved, e.g., by adding a color-neutral scalar-gluon 
coupling (pF^'^F^y to the QCD Lagrangian. The general singlet structure function F2 then reads 



F2,d) 



<e^>[ C{x, tts) ® q{x, tts) + Cy{x, a^)] (2.23) 



with 






C(f)q[x,as) C(j,g[x, as) 



1=0 



C,{x,as) = (?.''\l''''\] =f:a,maiC,M^). (2.24) 



G(fyy[X, CLs) 



1=0 



Here < e^ > is the average squared charge of the light quark flavors, and Cq the unit matrix 
times 6{1 — x). Using the scale dependence ( p.21| ) of the partons and Eq. ( |2.8| ) for the running 
coupling, the scaling violations of F2 can be written as 

"^^^ -- <e'>\c®k + pc'-{c®p® c-^ + I3C' ® c-^)® c; 



cilng2 ^^ ^ ,_ ^.. , ,^ T, V ^ ^ -. H yy ; ^ 7J 

+ {C®P® C-^+pC® C-^)® F2, (2.25) 

where the prime denotes the derivative with respect to a^. Both d F2/d\nQ'^ and F2 represent 
observables, hence the following combinations of splitting functions and coefficient functions are 
factorization scheme invariant: 

i"hom. = C®P® C-^+f3C'® C-\ (2.26) 

-finhom = C®k + l3C^- Jhom. ® C -, . (2.27) 

Putting back the perturbative expansions of all quantities involved, these invariants are given 



order- by-order in a^ [Bq, 34 , p3] by 



/horn. = a,Po + af{Pi+ Ci®Po--Po® Ci-/3oCi} + ... , (2.28) 

i"inhom = Ctcmfco + acmCtslfcl + C ik^ - Pq C^,i} + . . . . (2.29) 

The generalization to higher orders is straightforward if cumbersome. From these relations the 
changes AP and Ak of the splitting functions induced by a modification of the coefficients A C 



and A C^ can be easily determined. Recall, however, that a particular choice for the physical 
upper-row (electromagnetic) quantities C-^g, C^g and C^^ does not fully fix the transformation, 
as well-known for the hadronic DIS factorization scheme |^. In fact, we will use the DIS scheme 
for our subsequent numerical calculations in both the hadronic and the photonic cases. Starting 
from the usual MS scheme of fixed-order calculations [36-38,29], the transformation reads 



ACi 



r'Ms 



x] 



X) 



,^MS 




AC 



7,1 







(2.30) 



The lower-row choice in the hadronic part is the conventional continuation of the sum-rule 
constraint to all A^, that one of the photonic part is taken over from the DIS-y scheme of ref . [^ , 
see also Sect. 3.3. This concludes the general all-order framework, and we can now turn to the 
resummed anomalous dimensions and coefficient functions. 



3 Small-x resummations for the anomalous dimensions 
and coefficient functions 

In this section we briefly summarize relations for the resummed singlet anomalous dimensions 
in Mellin-A^ space which are used in the numerical analysis below. The anomalous dimension 
matrix is related to the corresponding singlet splitting functions by 



7(iV,a,) = -2P(iV,a,) 



dxx^ P(x, a. 



(3.1) 



The general form of the small-x resummed, unpolarized anomalous dimension matrix reads 

fe+i 



fc=0 ^^^ ^ m=0 



(3.2) 



with Us = CAOtsiQ'^)/^^- In the discussion below we include the LO [^ and NLO ^^ anomalous 
dimensions 7o ^^"^ 7i completely, and account for the small- a; resummed series in the leading 
(Lx, m = 0) and next-to-leading (NLx, m = 1) small- a; approximations 



k=2 \^^ ^ 



fc+1 



JO) 



(1) 



iT + N l\'^ + 0{N-) 



(3.3) 



3.1 The leading series 

The all-order resummation of the Lx series was performed in refs. 



lL{N,as 



"^Z^ I (0) (0) 

fc=0 \ gk,gq 9k,gg 



N 



-2 







Ct 



/Ca 1 )7L(iV,«.) 



with 7l(A^, as) being the solution of 



N 



P= — = 2V'(1) - ^(7l) - ^"(1 - 1l) = Xo{1l) ■ 
a. 



(3.4) 



(3.5) 



ip{z) denotes the logarithmic derivative of Euler's F-function. jl is a muhi-valued function for 



complex values of A^. The perturbative branch of the solution is selected by the requirement p9 



jL{N,as 






for lA^I ^ oo . 



(3.6) 



For small values of 7l(A^, ««) the asymptotic representation 



a.c 



7l = ^|1 + 2EC(2/+1)7l 



2/+1 



(3.7) 



,(0) 



holds, from which the coefficients g}. in Eq. (|3.4| ) can be determined iteratively. Here ({n) 
denotes the Riemann (^-function. Note that only (^-functions of odd integers contribute, which 
is expected for physical quantities in four dimensions, cf. ref 
were given up to k 



40| . Analytic expressions for g^ 



(0) 



14 in ref. 41 



99 



Later both analytic representations and the numerical 
values of these coefficients were determined to large values of k by various authors (cf., e.g. 

can be used as representation of 'Jl^N, a^ 



refs. [|^, ^). In fact, the series (|3.4|) can be used as representation of 7l(A^, a^) along typical 
integration contours (cf. Fig. 4 in Sect. 5) for the inverse Mellin transform, as shown in ref. [J3 
20 terms (or less) in Eq 



]4D are sufficient to obtain an accuracy of better than 0.01%. 
For the later numerical analysis, i.e., to perform the Mellin inversion of the solution of the 
evolution equations, it is necessary to locate the singularities of 7l(A^, a^) in the complex N 
plane. The singularities of the resummed leading singular part of the anomalous dimension 71, 
can be determined by differentiating eq. (|3.5|) with respect to p. 



The condition 



dp 

djL 
dp 



[-^'(7l) + ^'(1-7l)] 



^p'{-fl 



IT 



1 



2 sin (iT-fL) 







(3.8) 



(3.9) 



yields the value of the resummed anomalous dimension 7^, at the branch points. The correspond- 
ing values of p are then determined by Eq. ( p.5|) . For the perturbative branch one obtains 



7i - 
72,3 



1 
2 

-0.425214 ± 0.473898 i 



pi = 4 In 2 , 

P2,3 ^ -1.41048 ±1.97212 t 



(3.10) 



cf. refs. [H, P 
For Rep 



The behavior of the real and imaginary part of 7l(p), is illustrated in Fig. 1. 
4 In 2, the first branch point. Re 7l(p) forms a 'roof at 7^ = 1/2 for Imp = 0, 
which remains stable over some distance in Re p. The imaginary part becomes discontinuous. 
At smaller values of Re p. Re 7l(p) develops two symmetric minima, and for even smaller values 
two additional maxima. Both extrema finally form the two other branch points, Eq. ( |3.10 ). In 
Im 7i these branch points manifest as single extrema of the corresponding curve for Re p = const . 



3.2 Next-order corrections 

The coefficients 7^^^^ and 7^^^^ of the NLx series in Eq. ( p.3| ) were calculated in ref. ^ . Recently 
also the first terms for Tgg\k have been determined H, |10|- AH these quantities are analytic 



functions of 7^, hence they do not introduce new singularities. The energy-scale dependent 



„(i) 



(1) 



contributions to Tgg^k have still to be derived, and the terms Tgq^k are also unknown so far in 



7nl(^>«. 






7t" 



'^^r 2 m 1 \ 

Igq.^h 7sg,NL / 



(3.11) 



DIS 



In the DIS factorization scheme, the function 7gg,NL(A^, ««) given by 
_ 2as ^ (1) (ois 



g, 2 + 37L - 37! [g(l - 7L, 1 + 7l)]' 
677 3-27i 5(2 + 27L,2-27i) 



^(7l) 



fc=0 



7 



(3.12) 



with 



i?(7) 



r(i-7)xo(7) 

r(l + 7){-7X'o(7)} 



1/2 



exp 



7V'(1) 



^ Xo(C) 



A:=0 



at; 



(3.13) 



Here B{x,y) denotes Euler's Beta-function, and 7j^°nl represents the anomalous dimensions in 
the Qo-scheme |^, in which the factor R{jl) does not appear. 

The presently available contributions to the splitting function xPqg{x) in the DIS scheme and 
their convolutions with a typical gluon shape are shown in Fig. 2 for as = 0.2, i.e., Q^ ~ 20GeV^. 
The LO splitting function vanishes like x for x ^ 0, their NLO counterpart is constant for x ^ 0. 
The strongly rising NLx result [|| therefore dominates below x ~ 10~^. This dominance persists 
after the convolution below x ~ 10^'^, although here the differences are considerably smaller 
than for the splitting functions themselves. 

The contributions oc Nf [9,46-48] and the energy-scale independent terms oc Ca fTO], ^ of 
lgg,NL{N,as), cf. also 1^, have been calculated recently. As shown in ref. 0, 7gg,NL(A^, «s) can 



be obtained from the larger eigenvalue 7+ of the resummed anomalous dimension matrix 



1 1 I ^ 

7± = 2 (^11 + ^9) ± 2^(^99 - 7gg) + 47qg7OT, 



(3.14) 



where 



7+ 



7- 



733 + 
Iqq ~ 



7qg7gg 

735 - 739 

733793 
733 ~ Iqq 



I99 



c 



+ 7^793 + O {a1f{a, 



C 
C 



Im - 7^793 + C) [a^J{as/N) 



(3.15) 



These relations result from the fact that the quarkonic (upper-row) entries in Eq. ( p.4| ) vanish, 
unlike the gluonic ones. Furthermore one has 7_ ^ —{Cp/Ca) ■ (2as)/(37r) -Tp due to Eq. (|3.11| ). 
In the Qo scheme, the present contributions to 7+ are determined as the solution of 0, p!0[] 



a 



^[Xo(7+) + «sXi(7+)] 



which yields 



7I" 



Xi(7l) 



(3.16) 



(3.17) 



after a perturbative expansion. Finally Xi{l) reads 

Xi(7)=Xp(7) + Xp^(7)+xf(7), 
with 



(3.18) 



go.a 

asXi 



an, na, 



OisXl 



Nfas Cb 



TT 



cos(7r7) 2 + 37(1 — 7) 



TT Ca V sin(7r7) / 3-27 (1 - 27)(1 + 27) 



NfUs 
47r 



^ (xo(7) + Xo(7)) - 0X0(7) - 



vr 



3cos(7r7) 2 + 37(1 — 7) 



4 (X^(7) + X[,(7)) + f f - t) ^°^^) + (^^(3) + 



sin(7r7)y 2(1 - 27) (1 + 27)(3 - 27) 

Hi) 



TT 



37(1-7) 



TT \ cos(7r7) 
.sin(7r7)j 3(1 - 27) 



11 + 



7(1-7) 



;i + 27)(3-27), 
,(1) 



(3.19) 



The function /i(7) in eq. (p.l9| ), which contributes to ik,gg only for /c > 2, is given by 



M.)-|:".(ji^ + TT}T7) (3.20) 

in approximate form, with ai = 0.72, 02 = 0.28 and a^ = 0.16 [|T^. From these results the 
Qo-scheme anomalous dimension is then inferred by H] 



o,(i).Qo 



/3o 2 '^ 



r— = 7^'^ - — a'^^ln 7^°V-Y'n(7ia^^ 



7- 



(1) 



CV (1) ^ «sXi(7ff) 

n ,'19 ~ / / (0) N 



(3.21) 



W9 



Using the transformation of 7+ to the DIS scheme |^ , and employing Eqs. ( |3.12| ) and ( p.l5| ) as 
well, one finally arrives at 



7,,,NL - 7,,,NL+4^«, + [1 it:(7LJj7,,, 



gg.NL 



a 



' E [^/^Sf ^ + aZT + Ag^' 



k,gg 



fc=0 



a.. 



A^y 



The first terms for 75?mt" read : 



'99,NL 



7gg,NL 






+ 



23 as 

IH ^ + 

6 N 



71_7r2 
18 ~¥ 






— --7r2-8C(3)^' 
27 36 ^^ 'Ca. 



1276 71 2^79 7 , 52 C^. 

81 108 3 ^ ^ 120 3 ^ ^Ca. 






8384 233 2 284 91 4 ^ o/^.,^ 4 ^ 
162" +^^^')-720" +2C(5)--C(3)vr 



243 



^C3) - 16C(5)) 1; 



n) 



(3.22) 






+ 



/ 45928 638 



V 729 
31 



243 



-vr 



|C(3)^' - 2C(3)= 



^.' + H5c(5) + HH2c(3) 

2160 3 ^^ ^ 27 ^^ ^ 



rTT 



3024 



26 



^'C(3) - — C(5) - ^^C(3) - 80C(3)2J ^ 



n) 



(3.23) 
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A similar expression can be derived for 7fL^NL- Because of the yet approximate expression for 
h{'y) and the missing energy-scale dependent terms we will, however, only list the numerical 
values of those still preliminary expansion coefficients in Table 1. 

In Fig. 3 the different approximations to the splitting function xPgg{x) are displayed. Here 
both the LO and the NLO terms are flat for a; — >■ 0, while the Lx contribution |^ causes a strong 
rise as X ^ 0. Also shown are the presently known NLx terms just discussed. The addition of 
the quarkonic (NLxg^) contribution reduces the resummation effect almost down to the flxed- 



order results. The energy-scale independent gluonic terms [jlOl have an even stronger impact, in 
fact, they turn xPgg negative already at x ~ 10~^ for a^ ~ 0.2. A similar, but milder pattern is 
observed for the convolution x{Pgg ® /) with a typical gluon shape which illustrates the Q^-slope 
of the gluon density induced by Pgg. Note that the energy-scale dependent contributions to the 
NLx terms in Pgg have still to be calculated. These terms or yet unknown higher-order (NNLx) 
contributions may change the present behavior of Pgg{x). 

The leading singular contributions to the gluonic and pure-singlet quarkonic coefficient func- 
tions for the longitudinal structure function were also determined in ref. [H], 



Cf 



cl 




IL, 1 + 1l)Y 



27i, 2 + 27i) 



«(.).|^..£.^(|)' 



(3.24) 



(3.25) 



The numerical values of the first 20 expansion coefficients gi^gg, gl^gg in the DIS and Qo schemes, 
rk,c^, and coefficients contributing to g^J^^ and ^gk^gg are listed in Table 1 for complete- 
ness. These coefficients were calculated using the MAPLE-package [Q. The numerical values 



of 5^^^ /(41n2)'' were tabulated before for the DIS |^ and MS schemes [^. With a low num- 



ber of digits the values of gk^qg , "^k and c^ were given in ref. [^ as well. Either the direct 



expressions ( p.4| ), ( |3.12| ) or relations based on the corresponding expansion coefficients have been 
previously used in numerical studies M, |23|, p3. 



3.3 Photon-parton splitting functions 

Finally we have to consider the small-x higher-order corrections to the inhomogeneous terms 



kq = Pqj and kg = Pg^ in the photonic evolution equations (|2.21|) . These quantities arise from a 
subset of the diagrams leading to the gluon-parton splitting functions, with the incoming gluon 
replaced by a photon. Purely gluonic graphs do obviously not belong to that subset, as the 
photon can couple to the hadronic system only through qq emission. Hence, by comparing to 
the hadronic results discussed above, one obtains the most singular (S) small-x terms as 



ki 



^em kqfl[N- 



^0) + 



oo 

E 

1=1 



^fvn^ 



us 



k 



E 



Ciprri (l„ 



US 



(3.26) 



I.e., kq (up to its scheme-independent constant LO term) is definitely beyond the current NLx 
approximation, whereas kg can receive contributions at this order. Likewise, the constant but 
scheme- dependent term C^ i is NLx, and all higher-order photonic coefficient functions are be- 
yond that approximation. A 1/A^-term is indeed present in the known NLO result kgi. This 
term, however, vanishes after transformation to the DIS or DIS^ schemes. In fact, this cancella- 
tion can always be achieved at all orders, as we will show now. 
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k 




^il (Qo) 


^il (DIS) 


rk 






1 

2 
3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 


1.00000 E+0 
0.00000 E+0 
0.00000 E+0 
2.40411 E+1 
0.00000 E+0 
2.07386 E+1 
1.73393 E+1 
2.01670 E+0 
3.98863 E+1 
1.68747 E+2 
6.99881 E+1 
6.61253 E+2 
1.94531 E+3 
1.71768 E+3 
1.06433 E+4 
2.55668 E+4 
3.67813 E+4 
1.71685 E+5 
3.75379 E+5 
7.36025 E+5 


1.00000 E+0 
2.16667E+0 
2.29951 E+0 
5.06561 E+0 
8.79145 E+0 
1.90521 E+1 
4.58482 E+1 
9.24159 E+1 
2.31063 E+2 
5.59958 E+2 
1.24822 E+3 
3.25381 E+3 
7.93653 E+3 
1.89275 E+4 
4.98520 E+4 
1.23011 E+5 
3.06504 E+5 
8.07771 E+5 
2.02210 E+6 
5.17873 E+6 


1.00000 E+0 
2.16667 E+0 
2.29951 E+0 
8.27109 E+0 
1.49249 E+1 
2.92268 E+1 
1.02812 E+2 
1.94887 E+2 
4.85100 E+2 
1.52444 E+3 
3.11451 E+3 
8.58375 E+3 
2.47571 E+4 
5.47435 E+4 
1.56195 E+5 
4.26980 E+5 
1.01111 E+6 
2.89398 E+6 
7.69042 E+6 
1.91919 E+7 


1.00000 E+0 
0.00000 E+0 
0.00000 E+0 
3.20549 E+0 

-8.11742 E-1 
4.56248 E+1 
3.27070 E+1 

-2.95476 E+1 
1.08183 E+2 
3.99588 E+2 
1.33228 E+2 
2. 10243 E+3 
5.51142 E+3 
5.30316 E+3 
3.85296 E+4 
8.49086 E+4 
1.40384 E+5 
6.94998 E+5 
1.44307 E+6 
3.22738 E+6 


1.00000 E+0 
-3.33333 E-1 

2. 13284 E+0 
2.27231 E+0 
4.34344 E-1 
2.02643 E+1 
2.30315 E+1 
3.46449 E+1 
2.65004 E+2 
3.30038 E+2 
8.50371 E+2 
3.90849 E+3 
5.67433 E+3 
1.77680 E+4 
6.21982 E+4 
1.07028 E+5 
3.51475 E+5 
1.05058 E+6 
2.10341 E+6 
6.80747 E+6 


k 


„m{o.) (r\ \ 
9k,gg [HO) 


9Z? (Qo) 


gZ^;^ (DIS) 


gS^ (DIS) 


^9Z, 



1 
2 
3 
4 
5 
6 
7 
8 
9 

10 
11 
12 
13 
14 
15 
16 
17 
18 
19 


-1.00000 E+0 

-3.83333 E+0 

-2.29951 E+0 

6.42072 E+0 

-2.59764 E+1 

5.75787 E+0 

1.21690 E+2 

-2.66365 E+2 

5.43807 E+2 

1.96852 E+3 

-2.04998 E+3 

1.49302 E+4 

3.33837 E+4 

9.19579 E+3 

3.35804 E+5 

6.26484 E+5 

9.72892 E+5 

7.05626 E+6 

1.29507 E+7 

3.18568 E+7 


0.00000 E+0 
0.00000 E+0 
0.00000 E+0 
-1.19004 E+2 
0.00000 E+0 
-3.42186 E+2 
-2.28879 E+3 
-6.98786 E+2 
-1.11881 E+4 
-4.10835 E+4 
-3.39345 E+4 
-2.75933 E+5 
-7.55104 E+5 
-1.10387E+6 
-6.12763 E+6 
-1.45966 E+7 
-3.01102 E+7 
-1.30018 E+8 
-2.96814 E+8 
-7.45406 E+8 


-1.00000 E+0 
-3.83333 E+0 

-2.29951 E+0 
-6.04506 E+0 
-2.81814 E+1 
-2.60988 E+1 
-9.43607 E+1 
-3.54981 E+2 
-4.27828 E+2 
-1.67366 E+3 
-5.21390 E+3 
-7.99079 E+3 
-3.05607 E+4 
-8.37332 E+4 
-1.57171 E+5 
-5.64262 E+5 
-1.43675 E+6 
-3.14592 E+6 
-1.05144 E+7 
-2.59548 E+7 


0.00000 E+0 
0.00000 E+0 
0.00000 E+0 
3.96679 E+1 

-5.35750 E+1 
3.42186 E+1 
4.40583 E+2 

-7.39527 E+2 
1.11801 E+3 
4.86665 E+3 

-9.10195 E+3 
2.40902 E+4 
5.32758 E+4 

-9.58437 E+4 
4.46747 E+5 
5.92510 E+5 

-6.85258 E+5 
7.71985 E+6 
7.22515 E+6 
2.95797 E+6 


-1.65000 E+1 
0.00000 E+0 
1.48980 E+1 
-2.25291 E+2 
2.42631 E+0 
-2.09481 E+2 
-2.78211 E+3 
-2.70970 E+2 
-7.53012 E+3 
-3.82351 E+4 
-1.56381 E+4 
-1.73122 E+5 
-5.68231 E+5 
-5.13392 E+5 
-3.52577 E+6 
-9.13527 E+6 
-1.36171 E+7 
-6.85495 E+7 
-1.58648 E+8 
-3.23625 E+8 



Table 1: The numerical expansion coefficients in Eqs. (^, ( gl2D , (|3l3|) , (|3]2|) and (|]2|). 
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We proceed in two steps. We start by a purely hadronic transformation (like that one from 
MS to DIS), where G— * C = C + AC, and C^ remains unchanged. Then Eq. ( |2.27| ) implies 

= AJinh. = C- Afc + AC-fc. (3.27) 

We consider only such transformation terms, which are motivated in contributions to the actual 
electromagnetic coefficient functions. These are, however, of NLx order, as the k themselves, 
and hence one has Afc = on the NLx level: purely hadronic scheme changes modify the 
photon-parton splitting functions only beyond the NLx approximation. 

The second step is a purely photonic transformation, C^ — >• Cj = C^ + 6C^, with C 
untouched. Here Eq. (|2.27| ) yields 

= A /m, = C ■ Afc + /5 ■ a C; - /hom. ■ A C^ . (3.28) 

Solving for the changes Ak^ up to the order agm^r therefore leads to 



«1 = 


-thorn. ,0 ■ A C^_i 


«2 = 


— Ci ■ Aki + i"hom.,0 ■ A C^^2 + -fhom.,1 " A C 




m— 1 m—1 


m 


~ / . Cl ■ Ak-in-l + > ^ /hom.,1 ■ A C^^rn-l+1 , 




/=1 1=0 



7,1 



(3.29) 



where only terms have been retained which can potentially contribute to AfcNLx; if an NLx con- 
tribution to A C^ occurs in the transformation. By choice of the (unphysical) lower component 
of ACy^m the NLx pieces of km, arising, for example, in an MS calculation, can hence succes- 
sively be eliminated, without disturbing the vanishing of the upper term. At NLO, e.g., the lower 
component of zero in Eq. ( |2.3CI| ), as chosen in the DIS^ scheme |]32| achieves this cancellation. 
Therefore, without any loss of generality, vanishing resummed photon-parton splitting functions 
can be assumed at the NLx level. 

4 Less singular contributions 

The terms in the splitting functions Pij, which are less singular by one (or more) powers of 
ln(l/x) as X — i> than the leading contributions discussed in the previous section, are presently 
unknown in almost all cases. Such subleading contributions, however, can potentially prove to 
be as important as the leading terms, as also noted in a similar context in ref. [^. The splitting 



functions and coefficient functions enter observable quantities always via Mellin convolutions 
with the parton distributions. Since the parton densities are steeply rising towards small-x, but 
(at least in the hadronic case) small at large x, the structure functions probe the behavior of 
splitting functions and coefficient functions at medium and large values of x as well. 

The unpolarized singlet splitting functions are constrained by energy-momentum conserva- 
tion, see Eq. ( ^.IT] ). Also in other cases, however, as for the polarized singlet and the non-singlet- 
-|- evolutions, where no conservation laws apply, less singular terms with sizeable coefficients exist, 
for example in NLO, see, e.g., ref. |Q. In order to evaluate the possible impact of such terms 



in higher-order splitting functions, their numerical coefficients need to be estimated. At present 
the almost only source of information are the fully known LO and NLO splitting functions. 
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The dominant and subdominant terms in the small-x 1/A^-expansion of the LO and NLO 
singlet anomalous dimensions are recalled in Eqs. ( [4.1| )-( ^73D . In accordance with the main part 
of our numerical studies in Sect. 6, the results are listed for four light quark flavors. The first 
terms in the LO case are given by 

7gg,LO = +10.8793 A^- 6.82222 A^2^0(A^^), 

7gs,LO = -10.6667+ 11.5556 A^- 13.1852 iV2 + 0(Ar3), 

10 6667 
7ot,lo = 'j7 — + 8.00000- 9.3333 A^+ 10.0000 A^2 + 0(A^^), 

24.0000 9 . -^x .X 

Igg^LO = 77 — + 27.3333 - 5.1883 A^ + 17.0395 A^^ + 0(A^^) . (4.1) 

The corresponding expansions of the NLO anomalous dimensions read, in the DIS scheme, 



123.259 
277.333 

w 

91.2593 

1^ 
245.333 

1^ 



7?5lo = ^7 — + 405.863 - 684.836 A^ + 1197.52Ar2 ^ q(^j^3, 

7S!nlo = - ^' '^^'^^ + 846.222 - 1706.18 A^ + 2622.76 A^^ + 0{N^^ 

7Slo = + ''^tr^ - 453.512 + 809.030 A^ - 1344.89 A^^ + 0{N^) 

7Slo = +^^^^^- 988.210 + 2093.25 A^- 3109.08 A^^ + 0(A^^). (4.2) 



One notices that the first subleading terms occur with a sign opposite to that of the dominant 
one. Their prefactors are of the same order, but in most cases the subleading coefficients are 
by a factor of about 2 to 4 larger. At leading order the quarkonic terms are not singular as 
A^ -^ 0. The gg— term even starts proportional to A^, as a consequence of fermion-number 
conservation, Eq. ( ^^.161 ). The alternating structure continues towards higher powers in A^ with 



a similar pattern for the coefficients as observed for the first and second term. Note that this 
behavior is not a special feature of the DIS scheme, but is observed to a similar extent also in 
other schemes. As an example we give the corresponding coefficients also for the MS-scheme : 



94.8148 

1^ 



7Slo = — — + 253.026- 337.185 A^ + 623.259 A^2 + 0(A^^^ 



— 213 333 

7oTnlo = r^+ 461.449 -889.687A^+ 1501.16 A^2 + 0(A^^), 

an Q-| ^Q 

7^ NLO = +^ 361.805 + 658.108 A^ - 1048.43 A^^ + 0{N^) , 

7Slo = +^^^r^- 790.928 + 1616.55A^- 2423.77 A^2 + 0(A^^). (4.3) 



N 



In the structure function evolution, the difference between Eqs. ( [4.2|) and (|4.3| ) is compensated 
by the corresponding small-x terms of the coefficient functions. 

In the small-x resummation case, even partial results for subdominant contributions are only 
available for the gluon-gluon splitting function so far. The irreducible NLx contributions to 7^^ 



10 1 exhibit very large coefficients if compared to the Lx series M, see Table 1 and the comparison 



in ref. [0]. The introduction of terms with prefactors up to two times larger than those of the 
leading contributions, therefore, should provide conservative, non-exaggerating estimates for 
the possible impact of subdominant corrections. The following modifications of the resummed 
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anomalous dimensions beyond two-loop order, T{N,as), have accordingly been studied within 



refs. pT], E3|, pi p5 



(A) 


r(iv,«,)- 


-.r(iv,a,)-r(i,«,) 


(B) 


r(iv,a,)- 


-^T{N,as){l-N) 


(C) 


r(iv,a,)- 


-^T{N,as){l-Ny 


(D) 


r(iv,a,)- 


-^T{N,as){l-2N + N^). 



(4.4) 



The impact of the prescriptions (C) and (D) on the resummed NLx contribution to the 



splitting function xPqg{x,as) is illustrated in Fig. 2. At x = 10 ^ those terms reduce xP^ 



qg 



by factors larger than three, indicating the potential importance of less singular contributions. 
Also displayed in Fig. 2 is the convolution of Pqg with a typical hadronic gluon shape. The 
enhancement of the importance of non-leading terms by the Mellin convolution discussed above 
is obvious from the comparison of the two plots. 



5 Solution of the evolution equations 

In this section we derive the solution of the singlet evolution equations presented above. For 
technical convenience the analysis is performed in Mellin- A^ space where the convolutions turn 
to simple products. Recall that a unique analytic continuation of the anomalous dimensions 
from the integer moments to complex A^ exists |^. Thus a coupled system of two ordinary 
differential equations has to be solved at fixed A^. The x-space results are then obtained by 
a contour integral around the singularities of the final moment solutions f{N) in the complex 
A^-plane, e.g., that shown in Fig. 4. Due to /*(A^) = f{N*) it yields [0] 



xf{x) 



1 /""^ 

- / dzlm[e''^x-^f{N = C)] 

IT Jo 



(5.1) 



where C = c + ze^"^. For all cases considered here, c ~ 1 and ~ 37r/4 provide an efficient and 
numerically stable inversion. The latter choice of > it/2 leads to a faster convergence of the 
integral ( [5.1| ) as z — > cxd, see also ref. |^8[. At small-x, for example, a numerical accuracy better 
than 10~^ is easily achieved for upper limits as low as z^^x — 5. 



5.1 The general hadronic solution 

It is convenient to recast the evolution equations in terms of the running coupling a^ = as(<5^)/47r 
as independent variable, by combining the Q^ evolution ( |2.14| ) of the hadronic parton densities 
q with Eq. (|2.8|) for a^. Sorting the resulting r.h.s. in powers of a^, one obtains 



q{a,,N) 



dqia,, N) asPo{N) + a^,Pi{N) + alP2{N) + ... 



da. 



-al Po - alPi- af P2 



PqGs 



Pi 



p,{N) + a4P,{N)-t^Po{N) 



Pi 



P2 



-K^'<'^'-|^'"^'-{(l)-|hw)+'" 



1 

a.. 



Roim + Y.a'RkiN) 



k=l 



q{as,N) 



qia„N) 

(5.2) 
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Here we have simplified the notation by introducing the recursive abbreviations 



Po 

Rk — 'TrPk — / 'TrRk—i 

/3o fe /5o 



(5.3) 
(5.4) 



for the sphtting function combinations entering this expansion. As in Eqs. ( |5.3|) and (|5.4| ), we 
will often suppress the explicit reference to the Mellin variable N below. 

The splitting function matrices Pk of different orders k do generally not commute, especially 
one has 

[Rk>iiN),RoiN)]^0. (5.5) 

This prevents, already at NLO, writing the solution of Eq. ( |5.2| ) in a closed exponential form. 
Instead we proceed by generalizing the NLO method of ref. ||3^ to all orders]^ in a^. The 
corresponding ansatz of a series expansion around the lowest order solution. 



,LO, 



as,N) 



ao 



-Ro{N) 



q{ao, N) = L{as, Oq, A^) q{ao, N) , 



reads 



q{as,N) 



Uias, N) L{as, ao, N) U-\ao, N) q(ao, N) 

OO OO -. 

1 + ^ aj UkiN)] X(a„ ao, AT) [l + ^ a^ Uk{N)] ' q{ao, N) . 



(5.6) 



(5.7) 



k=l 



k=l 



The third, a^-independent factor in Eq. ( |5.7|) has been introduced to normalize the evolution 
operator to the unit matrix at Qq, instead of to the LO result (|5.6|) at infinitely high Q^. Inserting 
this ansatz into the evolution equations ( |5.2| ) and sorting in powers of a^ anew, one arrives at a 
chain of commutation relations for the expansion coefficients Uk{N): 



[Ui,Ro] = i?i+ C/i, 
[t/2,i?o] = R2 + R1U1 + 2U2, 



(5.8) 



k-l 



[Uk,Ro] = Rk + Y.Rk-iUi + kUk = Rk + kUk. 



i=l 



These equations can be solved recursively by applying the eigenvalue decomposition of the 
LO splitting function matrix, completely analogous to the truncated NLO solution with only 
Ui [Q, see below. One writes 

i^Q = r_e_ +r+e+ , (5.9) 

where r_ (r+) stands for the smaller (larger) eigenvalue of i?o, 



r± 



1 

Wo 



pf^ + pfJ ± V (^- - ^-T + 4^°'^^ 



(0) 
91 



(5.10) 



The first three orders were treated in a very similar manner in ref. 
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Ro-r^I , (5.11) 



The matrices e± denote the corresponding projectors, 

1 

e± = 

r± — Tip 

with / being the 2x2 unit matrix. Hence the LO evolution operator (|5.6|) can be represented as 

X(a„ao,iV) = e_(iV)U + e+(iV) p . (5.12) 

Vflo/ Veto/ 

Inserting the identity 

Uk = e C/fce_ + e [7^6+ + e^ C/fce_ + e^ U^e^ (5.13) 

into the commutation relations ( |5.8| ), one finally obtains the expansion coefficients in Eq. ( ^.7| ): 



^' = -* 



Ir 



e_Ri,e + e+Rue, 



e+ii^e e_Rj.e 



r_ — k 



(5.14) 



This relation completes the general structure of the hadronic singlet evolution. Note that the 
poles in Uk{N) at A^-values where r-{N) —r+{N) ±k vanishes are canceled by the U~^ term in 
the solution ( p.7|) . We are now ready to consider the presently available fixed-order and small-x 
resummation approximations. 

5.2 Fixed-order evolution 



In fixed-order perturbative QCD the expansion ( |2.15D of the splitting functions in powers of 



the strong coupling a^ is truncated at some low order k. Practical small-x calculations are 
presently restricted to NLO {k = 1), as the NNLO splitting functions P2{N) are not yet known 
for arbitrary values of A^, unlike the 2-loop coefficient functions P7| , [55[] and the /3-function 
coefficient (32 fl^]- Hence we confine ourselves to the NLO evolution here, the generalization to 



higher fixed orders being obvious. I.e., we keep the full results up to fc = 1 and put in Eq. (| 

Pk>2iN) = 0, (5.15) 

I3k>2 = 0. (5.16) 

The coefficients Pk>2 are also removed, for only all three quantities Pk, Ck and jSk together form 
a scheme independent set for the evolution of physical quantities like the structure functions 
^2(0;, Q^) or, in the case of polarized scattering, gi^x.Q"^). 

Two natural approaches have been widely adopted for the solution of the resulting NLO 
evolution equations. First one can solve Eq. ( |5.2| ) as it stands after inserting Eqs. ( |5.15| ) and 
(|5.16| ). Then still all orders in a^ contribute there and in the solution (|5.7l ), with the only 
simplification that the splitting function combinations (^.4|) are now given by 



This procedure is equivalent to a simple iterative solution of the system ( |2.14| ) and ( p.8| ), trun- 
cated aX k = 1. That technique is widely used in parton density analyses, e.g., in refs. [60-62]. 
The second approach uses power counting in a^ at the level of the evolution equation ( |5.2| ). 
There the a1 term in the square bracket involves P2 and (32 and can thus by considered as beyond 

17 



the present approximation. Consequently only the constant and the linear terms in as are kept, 
instead of Eq. ( |5.17| ) leading to 



^ NLO' ^ ^ NLO = 1 (p, - |p„) , Rl^O' ^ _ (5,18) 

Po ^ Po ^ 

In this approach it is furthermore natural to truncate also the evolution matrix [/(a^) after 
the linear term, since P2 would enter the determination of t/2 in Eq- ( ^-81) as well. Recall 
also that the final multiplication with the NLO Wilson coefficients does only cancel the scheme 
dependence of the linear a^ term in the evolution of the structure functions ( |2.25| ). Finally one 
can expand also U~^{ao) to first order in Eq. (|5.7|) , although this is not necessary, resulting in 

qS'^0(a„ N) = [i(a„ ao, N) + a, t/i(iV) L(a„ ao, N) - a^ L{as, ao, N) U^iN)] q{a^, N) , (5.19) 



where C/i is given by Eq. ( [5.14| ) with Ri = Ri of Eq. ( |5.18| ). This is the well-known truncated 



analytical NLO solution ^3| which has been employed, for instance, in refs. 

These two approaches obviously differ in NNLO terms only. The former procedure introduces 
more scheme- dependent higher order terms into the evolution of structure functions like F2{x, Q"^) 
or gi{x,Q'^) in a general factorization scheme. On the other hand, the latter method does not 
solve the evolution equations (|2.14|) literally, but only in the sense of a power expansion, i.e., up 



to terms of order k > 2. Therefore the first approach may be considered more in the spirit of 
the parton model, whereas the second is closer to a manifestly scheme independent expansion 
for physical observables. 

5.3 Small-x resummed evolution 

The resummed evolution of the parton distributions includes, to all orders in a^, the most singular 
small-x contributions to the splitting functions Pk- This inclusion is performed in the orders 
beyond the known fixed-order results. Thus the complete expressions for Po{N) and Pi{N) 
are used also here, and the difference to the previous section is restricted to the higher-order 
matrices Pk>2- Our notation in this section will directly apply to the evolution of unpolarized 
quark and gluon densities. Most of the subsequent discussion can, however, be easily transferred 
to the polarized singlet evolution [^ by replacing 1/N^^^ by 1/{N +lY^^^ with correspondingly 
modified coefficient matrices in all expansions. 

In the present case the most singular small-x terms in the evolution equations ( p.l4| ) behave 



like 1 / X all'^^ Ivl^ X [Q and 1/x aJl'^^ln^ ^x [§] as discussed in Sect. 3. In Mellin-iV space these 
additional resummation contributions replacing Eq. ( p. 15 ) read 



p Lx p NLx 

n>2(Ar) = ^ + z^,^. (5.20) 

In particular, the matrix P]^^ is related to the expansion coefficients gj^^^g in Table 1 by 

P^=(4C.)-sH,(^^;^j). (5.21) 

znl indicates whether only these leading small-x pieces are taken into account {Lx resummation), 
or whether also the next terms in Eq. (|5.20|) are kept (NLx resummation), 

_ J for Lx resummation , . 

1 1 for NLx resummation . 
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Recall that only the upper row of the NLx matrix is completely known at present P] . Results 
including that part only will be marked by NLxq below. 

With respect to the solution of the evolution equations, the situation is analogous to the 
fixed-order case, with the expansion parameter a^ replaced by N at each order fc > 2 of the 
strong coupling. The first option is obviously again the direct solution of Eq. ( p.2| ), now after 
inserting Eqs. (|5.16| ) and ( ^.20|) . To elucidate the generalization of the NLOtr. procedure to 
the resummed evolution, consider the splitting functions contribution Rk>2 ( p-4D arising from 
Eq. ( ^:20D : 



RlUN) 



1 
+ 



pr + ^NLivp: 



NLx 






NLx 



1 









> Lx 

fc-2 



+ ... 



^NL iVPf^2^ 



(5.23) 



The omitted terms involving higher powers of f3i and /?2, or (3k>3, are obviously even less singular 
as the last line for N —* 0. Therefore, if the power-counting in A^ is done on the level of Eq. ( |5.2| ), 
one immediately arrives at 



Rium 



1 1 



Pk^' + tNL 



n(^p 



NLx 
k 






1 T^Lx 

fc-1 



(5.24) 



Note that in the NLx (Lx) case the /3-function coefficients Pk>2 {Pk>i) do not contribute any 
more, and that /3i occurs only linearly in the former case. Thus Pi, which does not exhibit an 
Lx contribution, does not enter Eq. ( |5.24| ) in the present unpolarized case. All this is completely 
analogous to the R matrices (|5.18|) for the NLOtr. evolution. 

Before we turn to the U matrix for this second procedure, it is instructive to consider a small- 
X approximation to the unpolarized Lx evolution in this approach. Unlike in any other QCD 
singlet case, including the polarized leading small-x resummation |^, the splitting function 
combinations Rk,k'>2 do commute here: [Rk{N),Rk'{N)] = 0, due to the simple structure of 



the matrix (|5.21|) . This is still not sufficient for a closed solution of the evolution equation 
unless one also keeps the leading small-x contributions to Po,i(A^) only. Then, however, one 
eigenvalue of Pq vanishes, resulting in 



-px- 
-^0 




4CU 



(5.25) 



Using Eqs. (|3.4|) one obtains with this additional approximation 



9LT™^-(as,iV) 



exp 



2kr^"^^^^'''^^ 



exp 



/12L~ 



\PoN, 



1+E 



di{as,ao) 



1=1 



N^ 



q{ao,N) 



(5.26) 



analogous to the complete resummed non-singlet solution, see refs. [^ . This simple approximate 
expression, however, does of course not lead to any quark evolution. Eq. (|5.26|) can be completely 



transformed to x-space, cf. ref. |43 



^approx.j^^^^^^ = -F(as,ao,a;) ® el~"° q{ao,x) 



(5.27) 



19 



with 

F(as,ao,x) 

where 



S{l-x) + 



12L 



(3o\og{l/x) f^^ 



'^oMl/x) 
12L 



(l-l)/2 



Ii-iiz) 



12L 



log 



.1/2 



log 



aoJ 



(5.28) 



(5.29) 



and Iu{z) denotes the Bessel functions of imaginary argument. Similar expressions, e.g. in the 
double-logarithmic approximations, have been studied in detail long ago |55[ and were also 



considered recently [|66|. As compared to the complete Lx solution, however, the approximation 
( |5.27| ) yields gluon densities which are typically too large by a factor of 2 for an evolution from 
4 to 100 GeV^. Therefore we will not apply this approach in the following. 

We now proceed with the general resummed solution corresponding to the truncated NLO 
treatment where Uk>2 = 0, see Eq. ( |5.19| ). The generalization to the present case is to keep only 
those terms of Uk>2, which arise from the Lx and NLx pieces of Ri and of Rk>2 in Eq- (|5.23| ). 
Hence the NLO coefficient U[ = Ui is supplemented by 

[ f/'2, Ro] = Rr' + RfUf + 2U'2 (5.30) 

etc. Here R^ denotes the small-x contribution of Ri 

1 1 

and the corresponding expansion coefficient U-^ is given by 

[UlRo] = Rf+ C/f. (5.32) 

The final step analogous to the NLOtr. method is to keep the non-(N)L2; parts of t/i only linearly 
also in the inverse matrix U~^ of Eq. ( |5.7| ). This leads to 

"" 1 



RfiN) 



Pl'^" + ^NLiV(Pl 



NLx 



Pip 



Lx 



(5.3i; 



U'-\a,N) = l + aUf{N) + Y,a''Ui{N) + a{Uf{N) - U[{N)} 



k=2 



U-s^{a,N)-{l-aUl{N)} + l-a Ui{N) . 



(5.33) 



The last two terms represent the truncated NLO contribution. Insertion of this decomposition 
into Eq. (|5.71 ) finally yields (with L = L{as, Oq, N) for brevity) 



+ 



q[r(a„iV) = [L + asUiiN)L-aoLUiiN)]q{ao,N) (5.34) 

Usias, N) L Us\ao, N) - L - a, C/f (iV) L + a,L C/f (iV)] q{a^, N) . 
The first line is the NLOtr. result ( ^.19| ), the second line represents the resummation correction. 



5.4 The photonic solution 

We now turn to the parton distributions of the photon. In terms of the running coupling, the 
corresponding inhomogeneous evolution equation (|2.21| ) reads 

dq^ias, N) a,^{ko{N) + a,fci(iV) + a^,k,{N) + . . .} 



da. 



-a] Po - ai Pi - at P2 ■ ■ ■ 



+ had. 



K,{N) + J2a^^Ki{N) 



+ had. = K{as,N)+ had. . (5.35) 
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Here a^ra = a/47r denotes the electromagnetic fine structure constant, and analogously to 
Eqs. ( |5.3| ) and ( |5.4| ) we have introduced 

Ko = ^ko, (5.36) 

Po 

Ki ^ ^h-J:^K,_,. (5.37) 

Po j=i Po 

Finally 'had.' stands for the r.h.s. of the hadronic evolution equation ( p.2| ), with q replaced by 
q^ . The homogeneous component, qhom.) of the solution of Eq. ( ^.35|) is as derived in Sect. 5.1 
- 5.3. Hence only the inhomogeneous part, qjnh. = (f — 9hom. with q-^^^iX'^Oi ^) = 0; needs to 
be discussed here. This solution can be represented in terms of the hadronic evolution operator 

W^) as 

q;^^{as,N) = C/(a„ AT) a7^"(^) r daa^^^^^ U-\a,N)K{a,N) . (5.38) 

For the iterative solutions the remaining integral can be performed numerically. In the 
truncated procedures, on the other hand, U~^{a,N) has been expanded in Eqs. (|5.19| ) and 
(|5.33| ). In these cases the NLO photonic splitting functions, k^ = PP , should be treated in the 



same way as their hadronic counterparts P/'' previously, reducing Eq. ( |5.37| ) to 

K; = Ki = 1 (fci - ^feo) , K;>2 = . (5.39) 

Po ^ Po ^ 

In the following we will confine ourselves to physical factorization schemes like the DIS^ scheme 

|32| or the DIS scheme, where the photonic coefficient function 6*2,7 ^^^ been absorbed into the 

quark distributions. In these schemes K{as,N) does not receive any Lx and NLx resummation 

corrections beyond the leading order, as discussed in Sect. 3. Thus Eq. (|5.39|) applies to the 

NLOtr. photon evolution as well as to the corresponding resummed case. 

Inserting U^l^Q,{as, N) = 1 — asUi{N) and the expansion ( |5.39|) into the inhomogeneous 



solution ( ^.38|) , the a^ integration becomes obvious and one arrives at [^ 



-^q^J^^' {a„ N) = -[l + a,Ui{N)]{l-^L)[l-RoiN)]-'Ko{N) 

-{1- L) R^\N) [K,{N) - U,iN)KoiN)] + Oia^) . (5.40) 

In the resummed case a numerical integration remains over the all-order part of U'~^ in 
Eq. (|;33D. Defining 

A,,,{a,,ao,N) = - r^L-\a,ao,N)(u-s\a,N)-l + aU',{N))Ko{N) (5.41) 

Jao a^ ^ ^ 

the solution is, again using L = L{as,ao,N), given by 



+ Us{as,N) 



1(1 -^L)[l- Ro{Nr' Ko{N) + {1- L) R,\N) C/f (iV) Ko{N) 
as ao 



a.. 



- [1 + a, Uf{N)] - (1 - ^ X) [1 - Ro{N)]-' Ko{N) - {1 - L) Ro\N) Uf{N) Ko(iV) . 
tts ao 

This relation completes the Mellin-iV solutions of the fixed-order and resummed, hadronic and 
photonic evolution equations. We are now prepared to investigate the quantitative impact of 
the various approximations, for both the splitting functions and the solutions, on the parton 
densities and structure functions. 
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6 Numerical results 

In the following we study the numerical consequences of the fixed-order and resummed evolution 
kernels on the evolution of structure functions and some aspects related to potential uncertainties. 
Despite the impressive amount of small- a; structure function data already collected at HERA 
|jl|, , the present investigation does not aim at a comparative data analysis. Such an effort would 
require quite some flexibility of the non-perturbative initial distributions, especially for the gluon 
density which is only rather indirectly constrained by measurements of F2 and F^. A detailed 
data analysis requires rather precise and independent constraints on the small-x behavior of 
the gluon density, which are not yet provided by current measurements at HERA. A thorough 
implementation of heavy flavor (charm) mass effects in the resummation framework would be 
required as well. These mass effects are non-negligible at small x, where the charm contribution 
to F2 and Fl is substantial, in spite of the very large hadronic invariant mass, W"^ S> 4m^, 
cf. refs. [^. Both of these issues lie beyond the scope of the present paper. Since some of the 



resummation corrections turn out to be very large one would like to know as well the next-order 
resummed corrections to perform a detailed data analysis. 

In the following, therefore, the impact of the various anomalous dimensions and Wilson 
coefficients is instead illustrated for fixed initial parton densities of the proton and the photon. 
Accordingly all calculations are performed using the same values for as{Q'^). Specifically, the 
NLO relation (|2:T0| ) is employed with Ajv^=4 = 0.23 GeV above Q'^ = ml = (1.5 GeV)^ and. 



by continuity of as(Q^), with AAr^=3 = 0.30 GeV below that scale. Above (below) Q"^ = m^ 
the evolution equations are solved for four (three) massless fiavors, respectively, with c{x, ml) = 
c{x, ml) = 0. The small effects of the bottom flavor are entirely neglected. All subsequent results 
are derived in the DIS scheme discussed above, with the truncated solutions of Sect. 5 chosen 
as default. Only the singlet resummations described in Sect. 3 are taken into account, since the 
non-singlet contributions are suppressed at small x in the present unpolarized case, see Fig. 5, 
and its resummation correction is very small, cf. Sect. 2.2. 

6.1 Proton structure: fixed-order evolution 

Let us flrst consider the leading and next-to-leading order evolution of hadronic parton densities, 
putting emphasis on the small-x region. As the value of A^ =4^ given above, the initial distri- 
butions for our proton studies are adopted from the MRS (A') global flt [^ at a reference scale 



Qo = 4 GeV^. For the present purpose, the most relevant feature of these input densities is their 
small-x behavior which has been constrained by previous HERA data: 

xgP{x,Ql) oc xJ:p{x,QI) oc x"°-^^ for x -> . (6.1) 

Recall that, unlike the gluon distribution, the DIS-scheme quark densities represent observables. 
The LO and NLO small-x evolution of xS^ and xg^ to Q^ = 10 and 100 GeV^ is shown in 
Fig. 5 together with the initial distributions. The LO curves have been calculated, as indicated 
above, using the NLO input densities and as values in Eq. (|5.6|). Hence they do not represent 
results of an independent leading-order analysis, but directly illustrate the importance of the 
NLO terms relative to the LO contribution in Eq. ( p.l9| ). One notices that the perturbative 
stability of the presently available flxed-order evolution is theoretically satisfactory also at very 
low values of x. For instance, the NLO/LO ratio amounts to 1.25 (0.87) for the singlet (gluon) 
density at Q^ = 100 GeV^ and x = 10~^. Furthermore the numerical differences between the 
expanded solution (|5.19|) and the iterative approach ( 5.17 ) to the NLO evolution equations can 
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be considered as absolute lower limits on the uncertainties due to the unknown higher-order 
splitting functions. These offsets reach 3% at x ~ 10^^, while amounting to less than 1% for 
X > 10~^, see also Fig. 6 belowf]. Thus one may roughly expect a 5 — 10% small-x effect from the 
3-loop anomalous dimensions, if fixed-order renormalization group improved perturbation theory 
remains the appropriate framework down to x>10~^. Such an estimate is also corroborated by 
studies of the factorization scale dependence of F2 at small x [pO[] . 

It is conceivable, however, that the NLO contributions to the small-x anomalous dimensions 
are untypically small (as, for instance, l/N"^ terms are absent in 71, cf. Eqs. (4.2) and Table 1). In 
this context it is instructive to study the convergence of (formal) small-x expansions of anomalous 
dimensions and Wilson coefficients into the series 



^(a„iV) = ^ 



1=1 






(6.2) 



already at the LO and NLO level, where the full results are available. Also shown in Fig. 5, 
therefore, is the NLxg approximation to the leading-order evolution, for which just the 1/A^ 
terms of 7^^ and jg^ are kept together with the leading A^^O constants in 7^° and 7^°, see 
Eqs. (4.1). Note that this procedure is close to the well-known double-logarithmic approximation, 
cf. Sect. 5.3. As can be seen from the figure, this first approximation is very poor: xS (xg) 
exceed the full LO results by factors of about 1.7 (2.2), respectively, rather uniformly in x at 
Q^ = 100 GeV^, without any appreciable sign of improvement for decreasing values of x. 

Hence the question arises how many terms in the small-x expansion are required for arriving 
at a reasonably accurate representation of the complete fixed-order results. Accordingly, Fig. 6 
displays the ratios S'^pp^'^VE^^" and ^approx./^fuu f^^ ^^le LO and NLO evolutions with an in- 



creasing number of terms taken into account in the expansion (|6.2|) of 70 and 71 (in NLO the 
complete expression for 70 has been employed for all curves). One finds that in general three 
to four non-trivial small-x terms, i.e., contributions up to N^Lx at LO and N^Lx at NLO, are 
needed to achieve an accuracy of better than 10%. The NLO situation is not a peculiarity of 
the DIS scheme chosen here, as a corresponding MS analysis using Eqs. (4.3) yields similar 
results. Note that an interesting pattern emerges in both fixed-order cases: the approximate 
results alternate around the exact values with decreasing amplitude. If such a pattern were to 
persist to higher orders in a^, a first reliable estimate of their possible impact could be derived 
once two more non-trivial terms in all small-x expansions were known. This aspect may be of 
relevance for the resummed evolution addressed in the following. 

6.2 Resummed evolution 

We now turn to the effects and the relative importance of the Lx [^ and NLx [8-10] higher- 
order contributions to the splitting functions discussed in Sect. 3. In the present subsection the 
momentum sum rule (p. 17 ) is restored by prescription (A) of Eq. (4.4), i.e., P^^ and P^'' are 



supplemented by appropriate 6{1 — x) terms at all orders k >2. This procedure is the one with 
the least impact on the small-x results. Without any subtraction the sum rule would be violated 
by about 1% and 6% at Q^ = 100 GeV^ for the Lx and NLxg resummed evolutions, respectively, 
of our MRS (A') initial distributions. 

^In previous comparisons [gO| ^ deviations of up to 8% were found between these solutions. These large 
effects originated in the unfortunate choice of a traditional approximate NLO expression for as{Q'^), showing 



that the representation of the NLO solution of Eq. (2.8) tends to be relevant at the current accuracy of the data. 
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The resulting evolutions of the singlet quark and gluon densities are compared to the NLO 
distributions in Fig. 7 for Q^ = 10 and 100 GeV^. The relative importance of the available 
gluonic (lower row) anomalous dimensions is illustrated in Fig. 8(a). Consider first the effect of 
the Lx corrections [Q. These terms exert an appreciable influence on the gluon evolution, but 
much less on the quark densities in the kinematic region covered by the figure. At x = 10~^ 
and Q2 = iqo GeV^ e.g., ratios of c;Lx/^nlo ^ i^A and E^^/E^^o = 1.07 are obtained. This 
pattern obviously arises from the matrix structure of the Lx kernel ( |3.4| ); only at higher scales 
the quark effect fully approaches the gluon enhancement. 

The inclusion of the NLx, terms [§, i.e. the upper row entries in Eq. (|3.11|) , leads to only 
a small additional effect on xg{x,Q'^). The impact of these terms on xS(x, Q^) is, however, 
exceedingly large, as already evident from Fig. 2. These effects have been illustrated before, 
cf. refs. [0, ^], partly using different parametrizations for the input distributions at the starting 
scale Qo- The resulting enhancement with respect to the NLO evolution amounts to a factor 
of 2.8, for example, at x = 10~^ and Q^ = 100 GeV^. This huge correction is indeed entirely 
driven by the quarkonic anomalous dimensions, as also illustrated in Fig. 8(a): any 'reasonable' 
change of the gluonic splitting functions affects xT, by at most about 10%. Since only this 
one resummation contribution is known for the dominant upper-row quantities, a theory based 
estimate along the lines of Sect. 6.1 is not yet possible for the resummed quark distributions, and 
hence for the most important structure function, F2. We will therefore resort to the sum- rule 
prescriptions of Sect. 4 in the next subsection. 

In the gluonic sector, on the other hand, the theoretical situation has been improved recently 



by refs. [^, |T^, see Sect. 3. We remind the reader that the latter findings for 'jgg, although 
indicative, are not final yet, since the so-called energy-scale dependent NLx terms have still 
to be calculated. The effects of the known next-to-leading contributions are also presented in 
Figs. 7 and 8(a). The well-established qq contribution to 7^^^ [^], which is not expected to yield 
the largest subdominant terms, already removes more than half of the Lx effects on the gluon 
density at x<10~'^ for Q^ = 100 GeV^, see Fig. 8(a). The energy-scale independent gluonic 
contribution |T^ overcompensates the enhancement by the Lx and NLxg terms at all x in Fig. 7. 



As expected from Fig. 3, these terms are that large that they even cause a sign change in the 
slope of the gluon evolution for x < 10"'^. It seems natural to expect that yet missing terms either 
in NLx or unknown terms emerging in higher orders correct this behavior again. Thus a first 
uncertainty band of the possible resummation effects on xg{x, Q^) seems close to completion. In 
this context it should be recalled that the NLx anomalous dimension matrix is not yet complete, 
as 7^^^ still remains uncalculated. Note, however, that 7^°-' has an impact of less than 10% on 
both xS and xg in the Lx and NLx^ evolutions, cf. Fig. 8(a). Hence 7^^-' is presumably not 
a major source of uncertainty, as one may expect a rather moderate effect of this quantity as 
welig. 

As in the NLO case of Sect. 6.1, the differences between the iterative and the truncated 
solutions, Eqs. ( |5.23| ) and ( |5.34| ), of the evolution equations should yield a lower limit on the 



uncertainty due to missing terms in the anomalous dimensions. In fact, if the present small-x 
resummations collected the most relevant higher-order terms, a reduction of these offsets should 
take place with respect to the NLO evolution. The corresponding results are depicted in Fig. 8(b). 
While staying on the same level as in the NLO case for the Lx evolution, the offsets increase 
significantly as soon as the NLx terms are included, in particular for the singlet quark density: 



^T^Viio d-vTncir>+ cii-Tr^r-t io alc:r^ omi-iT-inT-f ini^ Ktt- 4-Vi^ +o/-*^- fViof ^A ' 



This expectation is also supported by the fact that 7^^ does not contribute to the eigenvalues of the resummed 



anomalous dimension matrix up to the NLa: level, see Eq. ( 3.15 ). 
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ratios S'*'^'' ^y'^trunc d q£ ^p ^q about 10% are found. This decreased stability may point to a 
larger uncertainty of the huge NLx quark enhancement. 

6.3 Structure functions and less singular terms 

Less singular (subleading) contributions to the anomalous dimensions, i.e., terms which do not 
exhibit the leading A^— >0 behavior, are vitally important for the LO and NLO evolution at small 
X, as demonstrated in Fig. 6: three to four terms in the expansion (|6.2|) are required for a good 
representation. In higher orders of a^ the leading N —>-0 poles become more singular, but so do 
the subleading contributions, and the number of singular pieces increases. There is, therefore, 
no obvious reason to expect terms less singular in A^ to be unimportant at low x in all-order 
approaches. At the present stage of the theoretical development, however, one has to rely on 
reasonable estimates for obtaining a first impression of their possible impact. For this purpose, 
we employ the momentum sum-rule prescriptions (C) and (D) of Eq. (4.4) for all anomalous 
dimensions with only one all-order term known presently, 

7g2W~^i:U^){l-2N + N'^) for tj = qq,qg,gq. (6.3) 

Here a = 2[3] for the prescriptions (C)[(D)]. In 7^^ we adopt the presently known NLx contri- 
butions, which are taken from Table 1. Hence only the A^^ or A^^ terms in this quantity are 
adjusted according to Eq. ( p.l7|) . In view of the structure of the A^-expansions of the LO and 
NLO terms estimates like Eq. ( |6.3| ) are conservative, i.e., they might underestimate the present 
uncertainties. 

In Figure 9 the resulting singlet quark and gluon densities are compared at NLx accuracy to 
distributions evolved with prescription (A). The subleading terms of the ansatz (D) are sufficient 
to overcompensate the huge leading resummation effect on xS(x, Q^) slightly. E.g., the NLx*^^^ 
result falls about 10% short of the NLO distribution at x = 10"^ and Q^ = 100 GeV^ Note 
that even the difference between the prescriptions (C) and (D), arising from the replacement of 
parametrically small N^Lx by N^Lx terms in the quarkonic anomalous dimensions, proves rather 
appreciable. This situation is similar for xg{xQ'^), where the effects of the sum-rule induced terms 
are positive because of the very large negative 7^^-' entries, cf. Table 1. The order of the curves 
is different here as compared to xS, since 7^3 differs between the cases (A), (C), and (D) only 
in the third term of the small-x expansion, unlike 7^^ which dominates the quark evolution. 
Although definite conclusions cannot be drawn from these prescription-dependent results, they 



nevertheless indicate clearly that the 1/A^ expansion ( |6.2| ) behaves similar as in the fixed-order 
cases. 

We now turn to the proton structure functions F2 and F^. Their small-x behavior, as obtained 
from the parton densities just discussed, is displayed in Figure 10. Since our calculations are 
performed in the DIS scheme, F2 is very closely related to the quark singlet distribution at small 
X. Thus the left side of Fig. 10 exhibits a pattern very similar to Fig. 9. The longitudinal 
structure function, on the other hand, in addition involves the resummed coefficient functions 



3 ^ \e, nAry 



7^ \ A: 



4 /?> \ "^ 



Cl^ and Cl,i represent the leading and next-to-leading order [|^, ^ coefficient functions, and 
the gluon and pure singlet resummation coefficients c^ are given in Table 1 and Eq. ( |3.25| ). 
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The additional corrections due to the coefficients c^ are, in fact, very large at the lower Q^ 
values shown: even the cross-section positivity constraint F^ < F2 is violated for x<3 • 10~^ 
at Q2 ^ 4 GeV^ for the MRS (A') initial distributions. At high Q^ > iqo GeV^ the effects of 
the coefficient functions and of the parton evolution become comparable due to the decrease 
of as in Ci,{as.,N^. The size of the low-Q^ effect shown in Fig. 10 (upper dotted and dashed- 
dotted curves), however, requires sizeable corrections by yet unknown higher-order terms in 
the small-x resummation of C^ or a large adjustment of the input gluon density. In fact, 
also the coefficient functions can be expected to receive relevant subleading corrections which 
are unknown at present. In order to derive a ffist estimate on their possible impact, the F^ 
calculations have been repeated with 

42--c42(l-2Ar). (6.5) 

The results of these calculations are shown in Fig. 10 (lower curves). This moderate correction 
term leads to an even drastic overcompensation of the leading resummation effect at low Q^ . 
This shows that for a more detailed understanding of the small-x behavior of F^lxjQ"^) the 
next-order small-x resummed corrections are required. On the other hand, direct measurements 
of Fl{x, Q"^) [| by the HERA experiments could help to constrain the size of missing terms in the 
coefficient functions. 

6.4 Photon structure at small-cc 

We now address, finally, the small-x evolution of the parton densities of the real photon. As 
outlined in Sect. 2, this evolution includes a specific inhomogeneous ('pointlike') piece in addition 
to the homogeneous ('hadronic') component. Whereas the latter behaves rather similar to the 
proton's parton distribution considered in the preceding subsections, the former is completely 
calculable in perturbation theory up to its dependence on the starting scale Q^. As discussed 
in Sect. 3 one may study the evolution of the photon structure function in a DIS scheme, where 
the inhomogeneous part does not involve any new resummed splitting functions at the present 
level of accuracy. It does, however, probe the resummed hadronic evolution matrix in a specific, 
different manner, cf. Sect. 5, and thus provides an additional laboratory for studying the possible 
effects of the small-x resummations. 

The reference scale Qq takes a somewhat different character in the photon case than in the 
pure hadronic evolution. It is a free parameter for the solution of the evolution equations, 
still, but only for certain choices of Qg can the separation between the homogeneous and in- 
homogeneous pieces approximately reflect the physical decomposition into the non-perturbative 
component, induced, e.g., by vector meson dominance (VMD), and a perturbative contribution. 
In fact this physical decomposition leads to Ql < 1 GeV^ J70|, |7l]], for a recent overview cf. 



ref. IjT^. We therefore choose Qq = 1 GeV^ in the following, unlike the proton case of Sect. 6.1- 
6.3. At this scale we adopt the NLO photonic parton distributions of GRV [0, as this is the 
only available NLO set with a HERA-like small-x rise of the hadronic component. The low-x 
behavior of these singlet and gluon densities is not given by a simple power law (6.1), however, 
but can approximately be written as 

S^(x,g^) ^ a;-o.22 ^ gi{^j:^Ql) ^ a;-o.i3...-o.22 f^^ 10-4<x< 10-2. (6.6) 



The measurements of Fl{x, Q^), ref. [^ , are 'indirect' and correlated with the F2 measurement. Their present 
experimental errors are still large. More precise results are expected from the data of the 1997 HERA run. 
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Here the effective rising power of xg'^ decreases witli decreasing x, cf. ref. ^4 . 

Tlie fixed-order evolution of xS''' and xg'^ is recalled in Fig. 11. As in the proton case, the 
LO solution has been calculated using the NLO (DIS scheme) initial distributions and the NLO 
values for ^^(Q^). The NLO/LO difference is slightly larger than in Fig. 5 due to the larger 
values of the coupling constant involved. Also shown in the figure is the NLO hadronic VMD 
contribution which is suppressed (in particular in the quark case) at large x, but dominant in the 
small- a; regime: it still amounts to about 80% of the full results for x < 10^^ at Q^ = 100 GeV^. 
Therefore one may expect a similar rise of F2 as observed for Fj* at HERA [Q. We shall consider 
now how the resummation corrections affect this picture. 

Fig. 12 present the effects of the various resummed small-x terms [4,8-10] on the evolution of 
the singlet quark and gluon distributions. The full results and the inhomogeneous contributions 
for Ql = 1 GeV^ are separately shown. These effects are considerably larger than those in 
Fig. 7: the NLx/NLO ratios reach factors of about 8 and 2 here for xT,'^ and xg'^, respectively, 
at X = 10~^ and Q^ = 100 GeV^. The dominant source of this greater enhancement are again 
the larger a^ values implied by the lower choice of Q^. The inhomogeneous components are still 
suppressed, although they are even more affected by the resummation corrections, as factors of 
up to 15 and 8 are found for the NLx/NLO ratios of Sinhom. and (^inhom., respectively Q. Recall 
that these latter results do not depend on any non-perturbative input distributions. Note also 
that xf^inhom. is much less affected by the 7^^^ corrections [|, |10[, since the main 'driving term' 



of the inhomogeneous solution is the purely quarkonic quantity k, cf. Eq. ( |5.42 ) 



The possible effects of less singular terms, using the same momentum sum-rule prescription 
as in the proton evolution, are illustrated in Fig. 13 for F2 and xg"^. The general pattern for 
the total results is analogous to the purely hadronic case of Figs. 9 and 10. The relevance of 
subleading corrections, however, is even more enhanced than that of the leading terms: F2 falls 
far below the NLO calculation for the ansatz (D), and the breakdown of the gluon evolution in 
NLx^"^-* already takes place at x<10~^. On the other hand, the less singular terms are much 
less effective in -F'2,inhom. at small x. In hadron-like cases their importance is magnified by the 
convolution with the (soft) parton densities, cf. Figs. 2 and 3. Here, however, the function 
x/cq oc x[1 + (1 — x)^], which plays the role of an 'input distribution', is very hard. Hence -F2^inhom. 
comes closer to a local probe of the small-x splitting functions than any inclusive hadronic 
quantity. 

Nevertheless the inhomogeneous part remains much smaller than than homogeneous piece of 
F2 for most scenarios of Fig. 13. It should be noted, however, that the latter contribution may 
be suppressed down to about the NLO results by a different choice of xg'^. It is conceivable, 
therefore, that -F2,inhom. is much more important in the resummed evolution than in the fixed- 
order case discussed above. Indeed, a small-x F2 considerably greater than about 1.2 times 
the VMD expectation could be considered as a signal for the presence of large resummation 
corrections in the quarkonic anomalous dimension. A measurement of F2 in the small-x region 
will, however, presumably only be possible with the 67 mode of a future e^e~ linear collider 0. 
Another theoretically cleaner, but experimentally also very difficult probe would be the structure 
of highly virtual photons, where the non-perturbative VMD part is suppressed and the calculable 
part becomes more important, cf. refs. [j73|. 



*In contrast, the NLO/LO ratio not shown in the figure is on a rather normal leveL 

^The possible kinematic coverage and necessary detector requirements have been studied in ref. [73| 
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7 Conclusions 

The effects of the resummation of the Lx and the known NLx small-x contributions to the ffavor- 
singlet anomalous dimensions and coefficient functions have been investigated, in a framework 
based on the renormahzation group equations, for the DIS structure functions F2 and F/ as well 
as for the photon structure function F2 ■ In this approach direct comparisons are possible with 
studies of the scaling violations of these structure functions based on LO and NLO fixed-order 
calculations. In order to allow for the most flexible comparison of different approximations to 
the all-order evolution equations, their general analytic moment-space solution has been derived. 
The largest small-x corrections to the quark densities and F2 are due to the resummed 
quarkonic NLx-corrections f^, whereas the effect of the gluonic terms is marginal here. For 
the gluon density, the Lx corrections 0] are moderately positive. Both the quarkonic and 
the energy-scale independent gluonic |]TU|] parts of the NLx gluon-gluon anomalous dimension. 



on the other hand, cause negative corrections which are that large that they overcompensate 
the Lx-terms. In fact, the latter terms lead to negative values for the total splitting function 
xPgg{x,as) for as = 0.2 and x < 10~^. This behavior probably signals the presence of other 
large positive contributions, either due to the energy-scale dependent NLx-terms or originating 
in terms of NNLx or even higher order. 

Contributions of NNLx order exhibiting a similar behavior can as well exist in the case of the 
quarkonic anomalous dimensions and the coefficient functions. This is suggested, for example, 
by the expansion of the fixed-order anomalous dimensions in powers of 1/A^ which leads to a 
good approximation only after three to four terms. Different ansatze for potential less singular 
terms have been studied numerically, showing that even the exceedingly large corrections due to 
the quarkonic NLx-corrections can easily be removed again. 

The longitudinal structure function F^ is in addition affected by the small-x contributions 
to the coefficient functions Cl |^. For lower values of Q^ the corrections become that large that 
the positivity constraint Fl < F2 can be violated for conventional input parton distributions. 
However, also this resummed coefficient function is very sensitive to subleading corrections. 

All these aspects show that also the next less singular terms need to be calculated, despite the 
enormous work that has been carried out so far to derive the resummed anomalous dimensions 
and coefficient functions [4,8-10,48,49], before firm conclusions on the small-x evolution of singlet 
structure functions can be drawn. Since contributions which are even less singular than these ones 
may, even then cause relevant corrections, it appears indispensable to compare the corresponding 
results to those of future complete fixed-order calculations. There the medium and large x terms 
are fully contained up to the respective order in a^. If extended to higher orders, in fact, also 
the RGE-improved fixed-order perturbation theory still seems to remain a viable candidate for 
the theoretical framework in the HERA regime. 

The small-x evolution of the real photon's parton structure has been analyzed in the DIS 
scheme. It has been shown that this scheme can be defined, without loss of generality, in such a 
manner that the photon-parton splitting functions do not receive any higher-order resummation 
corrections at NLx accuracy. Nevertheless the photon structure function F2 can provide an 
additional laboratory for studying the possible effect of small-x resummations, as the charac- 
teristic, calculable inhomogeneous solution of the evolution equations probes the low-x hadronic 
anomalous dimensions in a unique way: it comes closer to a local probe of the small- a; quarkonic 
splitting functions than any inclusive hadronic quantity. Unfortunately, this particularly interest- 
ing contribution is likely to be dominated by the hadron-like vector-meson-dominance part which 
behaves completely analogous to the photon structure and hence introduces the same uncertain- 
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ties and limitations due to the interplay of the anomalous dimensions and the non-perturbative 
initial distributions. 
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Figure captions 

Fig. 1 The real and imaginary parts of the perturbative branch of 7^ [Q as a function of 
p = N/ols- The the dash-dotted hnes are the contours through the singularities, Eq. (|3.1(]| ). 

Fig. 2 (a) The cumulative effect of the available contributions on the splitting function xPqg{x) 
for Us = 0.2. The fixed-order results are supplemented by the small-x resummed NLx 
corrections [§ beyond NLO. Also shown are the modifications induced by the subleading 
terms (C) and (D) of Sect. 4. (b) As in (a) but for the convolution of Pqg with a typical 
shape of an hadronic gluon density. 

Fig. 3 (a) The cumulative effect of the available terms for splitting function xPgg{x) for a^ = 
0.2. The fixed-order results are successively supplemented by the higher-order small-x 
resummed Lx correction [Q, the qq contribution to the NLx term 0, and the gluonic NLx 
energy-scale independent terms |T^. (b) As (a) but for the convolution of Pgg with a shape 
of an hadronic gluon density. 

Fig. 4 Integration contour in the complex A^-plane for the inverse Mellin transformation ( p.l| ) 
relative to the locations of the singularities of typical initial parton distributions (full 
circles), and those of the fixed-order (crosses) and resummed anomalous dimensions (open 
and closed diamonds for different values of Q"^). 

Fig. 5 The small-x evolution of the proton's flavor- singlet quark and gluon distributions in LO 
and NLO perturbative QCD. Two approaches to the solution of the evolution equations are 
compared in the NLO case, cf. Sect. 5. Also shown is the result of a small-x approximation 
(NLxq, see the text) of the LO splitting functions. The behavior of non-singlet quantities 
is illustrated by the total valence quark distribution x{u.u+dy). 



Fig. 6 The effects of the successive small-x approximations (|6.2|) to the LO anomalous dimen- 
sions 7o (left) and to the NLO corrections 71 (right) on the low-x evolution of xTP and 
xg'P^ cf. Eqs. (4.1) and (4.2). All results are displayed at Q^ = 100 GeV^ relative to the 
respective full calculations presented in Fig. 5. 

Fig. 7 The resummed small- a; evolution of the singlet quark and gluon densities as compared 
to the NLO results. The Lx and NLx [8-10] contributions are successively included, 
with the momentum sum rule implemented via prescription (A) of Sect. 4. The results for 
Q^ = 10 and 100 GeV^ have been multiphed by the factors indicated in the plots. 

Fig. 8 (a) The impact of the resummed gluonic splitting functions Pgg and Pgg of refs. |^, ||, ^ on 
xE^ and xg^ at Q^ = 100 GeV^. The NLxg results of Fig. 7 have been chosen as reference. 
(b) The offsets between the iterated and the truncated solutions of the evolution equations, 
see Sect. 5, for the NLO case and various resummation approximations. In all other figures 
the truncated solutions have been employed. 

Fig. 9 The possible effects of subleading corrections to the resummed anomalous dimensions, 
exemplified by the momentum sum-rule prescriptions (C) and (D) of Eq. ( |6.3|) , on the 
small-x evolution of the proton's parton densities. The results using the 5(1 — x) subtrac- 
tions (A) are as in Fig. 7. The MRS (A') [^ initial distributions (transformed to the DIS 
scheme) have been employed as in all other proton figures. 
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Fig. 10 The small-x behavior of the resummed structure functions F2 and Fl for the parton 
evolutions shown in the previous figure in comparison with the NLO results. The upper 
Fl curves include the resummed coefficient functions Cl f^, cf. Eq. (|6.4| ), the lower ones 
illustrate the possible impact of a sub leading contribution to C^. 

Fig. 11 The small-x evolution of the photon's singlet quark and gluon distributions in leading 
and next-to-leading order, starting from the NLO parametrization of [^ at Qq = 1 GeV^ 
as in all following figures. The hadronic (vector-meson-dominance induced) components 
are compared to the full results at NLO. 

Fig. 12 The resummed small-x evolution of xYT' and xg'^ to Q^ = 100 GeV^ as compared to 
the NLO results. The Lx ^ and NLx [8-10] contributions are successively included, with 
the momentum sum rule implemented via prescription (A) of Sect. 4. The effects on the 
photon-specific inhomogeneous solution are displayed separately. 

Fig. 13 The possible effects of subleading corrections to the resummed evolution kernels, exem- 
plified by the momentum sum- rule prescriptions (C) and (D) of Eq. ( |6.3| ), on the small-x 
evolution of the photon structure function F2 and the photon's gluon distribution. 
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